Magnetic field dependence of valley splitting in realistic Si/SiGe quantum wells 
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The authors investigate the magnetic field dependence of the energy splitting between low-lying 
valley states for electrons in a Si/SiGe quantum well tilted with respect to the crystallographic 
axis. The presence of atomic steps at the quantum well interface may explain the unexpected, 
strong suppression of the valley splitting observed in recent experiments. The authors find that the 
suppression is caused by an interference effect associated with multiple steps, and that the magnetic 
field dependence arises from the lateral confinement of the electronic wave function. Using numerical 
simulations, the authors clarify the role of step disorder, obtaining quantitative agreement with the 
experiments. 



Qubits in silicon are leading candidates for scalable 
quantum computing, owing to their favorable and well 
studied materials properties f^*^ Indeed, because of its 
prominence in the electronics industry, silicon may be the 
best understood semiconducting material. However, as 
devices continue to shrink in size, approaching the quan- 
tum regime, important questions arise. Unlike direct gap 
semiconductors, the conduction band structure in silicon 
possesses six symmetric minima or "valleys" that are not 
at the Brillouin zone center. Consequently, the minima 
are degenerate, and must be described by a valley index 
which competes with the spin index as a relevant quan- 
tum number in the qubit Hilbert spaced Therefore, to 
construct spin qubits in silicon, it is necessary to lift all 
valley degeneracy. 

A silicon quantum well grown on the [001] surface of 
strain-relaxed silicon-germanium is under tensile strain, 
causing the four lateral valleys to rise significantly in 
energyi^ At low temperatures, only the two low-lying val- 
leys are populated. The remaining two-fold degeneracy 
can be removed by the sharp confinement potential of 
the quantum well interface.'^ Theoretical estimates sug- 
gest that the resulting valley splitting can be of the order 
of 1 meV « 12 K,^ which is sufficiently large for quantum 
computing. However, recent experiments in SiCe^i^iSiii^ 
measure a valley splitting much smaller than the theo- 
retical prediction. There is currently no explanation for 
this discrepency.^^ Indeed, a prevalent theory^^ predicts 
an enhancement of the valley splitting in a magnetic field 
that is different from the experimental observations. 

In this letter, we describe a single-electron valley split- 
ting theory for silicon quantum wells grown on a vicinal 
substrate, building upon an initial suggestion by Ando.^^ 
Such miscuts are often incorporated into Si/SiGe het- 
erostructures to ensure uniform growth surfaces and to 
avoid step bunching. The resulting quantum well, ob- 
tained by conformal epitaxial deposition, is misaligned 
with respect to the crystallographic z axis, as shown in 
Fig. [TJ We now describe the effective mass theory and 
explain how the presence of interfacial atomic steps sup- 
presses valley splitting. 

For silicon strained in the [001] direction, the effective 
mass wavefunctionM can be written as a sum of contri- 
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where the terms inside brackets are Bloch functions. The 
phase angle (p is determined by the position of the quan- 
tum well interfaceJi^ The two orthogonal valley states 
correspond to a phase difference of Acj) = tt, although 
the absolute value of is unimportant in the present dis- 
cussion. It is crucial to note that while the valley minima 
occur along the z axis at ±kQZ, the quantum well normal 
is tilted away from z. For a slowly varying confinement 
potential, the two valleys have the same envelope func- 
tion, and are essentially independent. However, sharp 
variations in the potential cause the effective mass ap- 
proximation to break down. Examples include the cen- 
tral cell potential near a shallow donor and the sharp 
band offsets at the interface of a quantum well.^^ Leading 
order corrections to the effective mass theory^^ give the 
valley splitting 



(2) 



where the coupling potential Vy{r) decays a few 
Angstoms from the interface. Because this decay length 
is so small, compared to effective mass length scales, it 
can be represented as a (5 function: 



Vvir) = VyS{z - z,), 



(3) 



where Zi{x) is the position of the interface along the x 
axis. The coupling parameter Vy contains atomic scale 
information that must be obtained from ab initio theo- 
ries such as tight-binding theoryji^ or from experiments. 
The phase factor in Eq. ^ reflects the valley separation 
of 2fco, and plays a crucial role in the theory by introduc- 
ing interference effects, as we shall see. Note that large 
modulation doping fields usually confine the wave func- 
tion to one side of the quantum well, so only one interface 
potential has been included in Eq. 

We can now explain the suppression of the valley split- 
ting. In a tilted quantum well, the interface position Zi{x) 
describes the atomic steps, as shown in Fig. [TJ The S 
function in Eq. ^ reduces the integral to a sum over 
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FIG. 1: Quantum well step geometry, with crystalographic 
axes (a;, y, z) and rotated axes {x ,y , z ), where y — y . 



steps with phase angles differing by 2kob ~ O.SStt. Here, 
ko ~ 0.85(7r/26), and b ~ 1.358 A is the atomic step 
height. Thus, the phases from consecutive steps interfere 
almost fully destructively. For a delocalized electron, the 
wavefunction extends over an infinite number of steps, 
resulting in the complete suppression of Ey . In the pres- 
ence of a magnetic field, the electronic wave function is 
confined to a finite number of steps, leading to a non- 
vanishing valley splitting. 

We can estimate the valley splitting at low magnetic 
fields from Eq. ([2]), assuming the idealized (but unre- 
alistic) miscut geometry shown in Fig. [U with uniform 
steps of equal width. A high-field tight-binding theory 
was presented in Ref. [l^ for the same geometry. A 
good approximation for the wave function of an elec- 
tron in a perpendicular magnetic field and a tilted quan- 
tum well is given by the well-known solution for a flat 
quantum well,^® in the rotated basis {x',y',z') shown in 
Fig. [1] For the lowest {n = 0) eigenstate correspond- 
ing to the symmetric gauge, A = {—y',x',0)B/2, we 
obtaini^ F{r') = F,y{x' ,y')F,{z'), where F,y{x',y') = 
e-i^'^+v'^)/iil /^2ttI% and F^{z') is the subband enve- 
lope. Here, the electronic wave function is confined over 
the magnetic length scale Ib = ^/h/\eB\. For large 
Ib, we can ignore the discreteness of the atomic steps 
in Zi{x). Performing the integral in Eq. we obtain 
Ey ~ 2u„F|(0)e-2('=o'B^)', where 6* < 27r is the miscut 
angle. Thus, for uniform steps, the valley splitting is 
suppressed exponentially in i3, as confirmed by the nu- 
merical analysis described below. 

The previous results are strongly affected by even a 
little disorder. At low fields, disorder can enhance the 
valley splitting by many orders of magnitude, showing 
that the exponential suppression occurs because of a del- 
icate cancellation of the phase terms in Eq. To 
investigate the effects of disorder, we have performed 
numerical simulations of the valley splitting. We intro- 
duce adjustable parameters describing the amplitude of 
normally-distributed fluctuations, or wiggles, of the step 
edges. The fluctuation model incorporates alternating 
smooth and rough step edges, as consistent with the ex- 
perimental data.'^^ Typically in our simulations, the am- 
plitude of the rough step fluctuations is set at its charac- 



teristic "large" value s, corresponding to the average step 
separation. We also introduce a step-bunching parame- 
ter, which allows neighboring step edges to overlap. Such 
bunching is known to occur at silicon growth interfaces, 
particularly under the influence of strain.^'' A represen- 
tative step profile is shown in Fig. [IJa) . To make contact 
with the data of Ref. 1^, we consider a 2° miscut. 

For a given, randomly generated step profile, we eval- 
uate Eq. ([2]) numerically. The resulting valley split- 
ting now depends on the electron position, since disor- 
der breaks the translational symmetry. In the symmetric 
gauge, the unperturbed magnetic wave function can be 
centered, degenerately, at different points in space. By 
computing the valley splitting for wave functions cen- 
tered at each of these points, we obtain a valley splitting 
landscape, Ey{r'), including peaks and valleys. Electrons 
are attracted to the valley splitting peaks in order to min- 
imize their total energy. We observe that peaks always 
occur near broad step fluctuations, or "plateaus," and 
that bunched steps tend to produce the strongest valley 
splitting peaks. By tracking the dominant peak in a given 
cell, we can determine Ey(B) for a fixed step profile. In 
Eq. @ we use Vy ~ 1.2 x 10~^^eV • m, as obtained in 
Ref. llSl. which corresponds to a Sio.yGeo.a/Si/Sio.TGeo.a 
quantum well. 

Some typical simulation results for Ey{B) are shown 
in Fig. [2Ib). By setting the step fluctuations to their 
characteristic "large" value, as described above, we ob- 
tain reasonable agreement with the experimental data of 
Ref. [13 for any value of the step bunching parameter. By 
fine tuning the step-bunching parameter, we obtain the 
best quantitative agreement when the average bunching 
length is in the range of 10-20s. 

In the intermediate field range 0.3 < i? < 3 T, the 
experimental data points in Fig.[2]are strikingly linear 
The corresponding simulation results exhibit a slightly 
sublinear field dependence. More generally, we find that 
the precise shape of the Ey{B) curves depends on the par- 
ticular disorder realization and the fluctuation model. It 
is very likely that other disorder models, not considered 
here, could produce different curve shapes, including the 
linear behavior of the data. One idealized model giving 
a linear valley splitting is the "plateau" model, corre- 
sponding to wide, localized, flat plateaus surrounded by 
relatively uniform steps. The valley splitting contribu- 
tion from the ordered steps surrounding the plateau is 
strongly suppressed, due to the cancellation effects de- 
scribed above, so that the plateaus dominate the valley 
splitting. The scaling theory for an isolated plateau is 
obtained by introducing the concept of an "excess area" 
A relative to a perfectly uniform step configuration. For 
a small magnetic field, the electronic wave function is 
spread out over the large area 2ttI^. The resulting valley 
splitting in this model is 



Ey ~ VyF^{0)A/TTl 



(4) 



which is linear in B. We can obtain an experimental 
estimate for A from the data of Ref. [T^. Using theoretical 
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FIG. 2: Realistic valley splitting simulations, (a) A ran- 
domly generated step edge profile, with alternating smooth 
and rough step edges. The circle shows the wave function 
radius Ib for _B = 3 T, centered at the position of the valley 
splitting maximum, (b) Simulation results for the valley split- 
ting (solid lines). Experimental data points (circles) are taken 
from Ref. fiol The top three curves assume the same "large" 
amplitude of step fluctuations, with different step bunching 
parameters. From top to bottom: (i) strong bunching, (ii) 
weak bunching, (iii) no bunching, (iv) fluctuation amplitude 
reduced by a factor of 5. Note that the top curve corresponds 
to the disorder realization in (a). 



estimates for Vy and ^^(O);^ we find that A ~ 18s^, 
where s = 3.9 nm is the average step separation for a 
2° miscut. This result in general agreement with our 
simulations. For a round plateau, this corresponds to a 
diameter of about 5s. 

The plateau scaling model must break down for very 
small fields (large Ib), when the wave function encloses 
additional plateaus. At high fields, the scaling also breaks 
down when the wave function is confined to a single step, 
saturating at the theoretical upper bound of Ref. H. For 
a 2° miscut, this crossover occurs at about 25 T. Thus, 
the scaling expression in Eq. and possibly the linear 
experimental data, may correspond to crossover behav- 
ior. Finally, we note that the magnetic field dependence 
of the valley splitting should be first order in the valley 
coupling parameter u„, consistent with Eq. since it 
involves breaking the positional degeneracy of the mag- 



netic (Landau) eigenfunctions. 

Valley states can be detrimental for spin-based quan- 
tum computingi^ Large valley splittings are needed 
to avoid thermal excitations outside the spin-1/2 qubit 
Hilbert space associated with the valley ground state. ^ To 
accomplish this in an experimental setting, the present 
results suggest that we should use substrates without 
miscuts, although it is usually difficult to eliminate large 
wavelength roughness in conventional, strained devices. 
An alternative approach is to develop heterostructures 
utilizing strain-s/iarm;; techniques, which can, in prin- 
ciple, provide step-free interfaces. Indeed, large valley 
splittings have be obtained in similar unstrained silicon 
oxide structures.— 

Finally, we point out that a tilted magnetic field can be 
used to test the proposed mechanisms for valley splitting 
suppression. When the magnetic field is tilted away from 
the growth axis {z'), we would expect strong variations 
in the valley splitting, depending on the relative angle of 
the field with respect to the quantum well. 

In summary, we have demonstrated that experimen- 
tal observations of the magnetic field dependence of the 
valley splitting are consistent with a theory in which in- 
terference between different atomic steps at an interface 
causes the valley splitting to decrease, as the size of the 
electronic wave function increases. The present results 
suggest that valley splitting can be increased by using 
substrates without miscuts, although it is usually diffi- 
cult to eliminate large wavelength roughness in conven- 
tional, strained devices. Alternatively, valley splitting 
can be controlled by constraining the wave function with 
electrostatic gates, as demonstrated in Ref. [23 . 
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